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ABSTRACT 

By implementing widely-used equations of state (EOS) from Lattimer & Swesty (LS) and H. Shen 
et al. (SHEN) in core-collapse supernova simulations, we explore possible impacts of these EOS on 
the post-bounce dynamics prior to the onset of neutrino-driven explosions. Our spherically symmetric 
(ID) and axially symmetric (2D) models are based on neutrino radiation hydrodynamics including 
spectral transport, which is solved by the isotropic diffusion source approximation. We confirm that in 
ID simulations neutrino-driven explosions cannot be obtained for any of the employed EOS. Impacts of 
the EOS on the post-bounce hydrodynamics are more clearly visible in 2D simulations. In 2D models 
of a 15 Mq progenitor using the LS EOS, the stalled bounce shock expands to increasingly larger radii, 
which is not the case using the SHEN EOS. Keeping in mind that the omission of the energy drain by 
heavy-lepton neutrinos in the present scheme could facilitate explosions, we find that 2D models of an 
11.2 Mq progenitor produce neutrino-driven explosions for all the EOS under investigation. Models 
using the LS EOS are slightly more energetic compared to those with the SHEN EOS. The more 
efficient neutrino heating in the LS models coincides with a higher electron antineutrino luminosity 
and a larger mass that is enclosed within the gain region. The models based on the LS EOS also 
show a more vigorous and asphcrical downflow of accreting matter to the surface of the protoneutron 
star (PNS). The accretion pattern is essential for the production and strength of outgoing pressure 
waves, that can push in turn the shock to larger radii and provide more favorable conditions for the 
explosion. Based on our models we investigate several diagnostic indicators of the explosion that have 
been suggested in the literature, e.g., the amplitude of the standing-accretion shock instability mode, 
the mass weighted average entropy in the gain region, the PNS radius, the antcsonic condition, the 
ratio of advection and heating timescales, the neutrino heating efficiency, and the growth parameter 
of convection. 

Subject headings: equation of state — hydrodynamics — neutrinos — stars: neutron — supernovae: 
general 



1. INTRODUCTION 

Core-collapse supernova explosions are triggered by the 
gravitational energy released during the transition from 
a stellar core to a protoneutron star (PNS). The thermo- 
dynamic conditions obtained at the center of the stellar 
core are temperatures on the order of tens of MeV and 
densities on the order of normal nuclear matter density 
(3 x 10 14 g cm -3 ). There are only few equations of state 
(EOS) available for the se co nditions. Several st udies by 
iTakahara fe Satol (| 19821 ) and lBaron et~aT1 (|1985l ) were pi- 
oneering, investigating the impact of the EOS on the ex- 
plosion dynamics, albeit in a phenomen ological approach 
(see collective references in|Bcthc 1993). The most com- 
monly used nuclear EOS in recent supernova simula- 
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tions are the EOS from lLattimer fc Swestvl (|199lD (here- 
after LS) , based on the incompressible liquid-drop mode l 
including surface effects, and from IShen et all (Il99l 
(SHEN). The latter is based on relativistic mean field 
(RMF) theory and the Thomas-Fermi approximation. 
R ecently, a new nuclear RMF EOS has been published 
in lHempel fc Schaffner-Bielichl (|2010f ). which in addition 
can be obtained for several nuclear paramctrizations and 
provides detail ed information about the chemical comp o- 
sition (see also lShen et al.ll201lHFurusawa et al.H201lT ). 

Thus, extensive studies have been performed to in- 
vestigate the impact of the nuclear EOS on supernova 
simulations. But most of these simulation s are lim- 
ited to spherical symmetry. iThompson et al.l (|2003[ ) in- 
vestigated the effects of the incompressibility parame- 
ter on the evolution of temperature and neutrino lu- 
minosities during the postbounce evolution based on 
the three different incompressibilities of the LS EOS 
(A" = 180, 220, 375 Me V). Only a sma l l impact on 
the dynamics was found. ISumivoshi efaLl (l2006l . l2007h 
performed lo ng-term simulations of a 4 Mq progeni- 
tor star from IWooslev fc Weaver! (fl995h comparing LS 
(A" = 180 MeV) and SHEN, focusing on the evolution 
unti l black-hole formatio n and the emitted neutrino sig- 
nal. iFischer et alj (|2009[ ) included in addition progenitor 
stars from different stellar evolution groups in the mass 
range of 40-50 Mq and studied the neutrino signal and 
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collapse to a black hole in a syste matic way. Their results 
were confirmed and extended bv lO'Connor fc Ottl (|201lD 
applying a simplified treatment of neutrino transport. In 
iHempel et al.l (|2012t ). the authors applied their new EOS 
and systematically explored the hydrodynamic evolution 
and emission of neutrinos for different nuclea r parameters 
within the class of RMF EOS . Moreover, iSagert et al.l 
(|2009t ) and lFischer et al.l (|2011h explored the appearance 
of quark matter via a first order phase transition in su- 
pernova simulations during the early post-bounce phase, 
which triggers the explosion even in spherically symmet- 
ric models. 

Furthermore, EOS studies based on multi-dimensional 
supernova models have become availab l e rec ently. 
iKotake et ail (|2004D and iScheidegger et al.l (j2010l) per- 
formed two-dimensional (2D) and three-dimensional 
(3D) simulations, respectively, and investigated the ini- 
tial post-bounce phase and gravitational wave emis- 
sion focusing on effects of rotation and convection with 
simplified neutrino treatments. They found that the 
gravit ational wave signal depends strongly on the EOS 
(e.g.. | Kotakei 120111 for a recent review). In addition, 
iMarek eta l. (2009) investigated the EOS dependence of 
high frequency variations of the neutrino luminosity and 
gravitational wave emission by performing 2D simula- 
tions that include detailed neutrino-radiative transfer. 

Several 2D studies obtain neutrino-driv en explosions 
of massive iron-core progenitors (see e.g. IMarek et al. | 



20091: 



Bruenn et al.ll2009t ISuwa et al.ll2oTol : IMiiller et alJ 



2012al) . aided by the standing accretion-shock instability 



(SASI) and neutrino-driven convection. These simula- 
tions were evolved until the expanding explosion shock 
left the central core (^1000 km) between about 400- 
600 ms after bounce (depending on the progenitor and 
details of input physics) . Most multi-dimensional super- 
nova studies used the LS EOS with the inco mpressibility 
of K = 180 MeV, which is a very soft E OS dBuras et al. 
I2006al: ISuwa et al.ll2010D . In addition, IMarek fc Janka 
(120091) performed 2D s imula tion using the stiffcr EOS 



from rHiilebrandt et alJ ()1984l ) , and found that the softer 
EOS results in a more optimistic situation for a possible 
onset of an explosion. However, their simulation with 
the stiffer EOS had been performed for a shorter post- 
bounce time than the softer EOS. Hence, a final con- 
clusion cannot be drawn@ Even further, they employed 
only two EOS, each of which has a different incompress- 
ibility and symmetry energy. Note that the entire set of 
nuclear parameters, not only the incompressibility and 
the symmetry energy, but also their density dependence 
determine the resulting EOS. 

In this paper we present results of numerical simu- 
lations of core-collapse supcrnovae of massive iron-core 
progenitors. The spectral neutrino transport is treated 
by the Isotropic Diffusion Source Approximation (IDSA) 
(ILiebendorfer et al.H2009h . We employ four EOS of LS 
(with the three incompressibilities) and SHEN. We ap- 
ply these different EOS in axially symmetric (2D) sim- 
ulations and investigate the differences obtained. In the 
analysis of the data we focus on the post-bounce phase 

8 Employing a stiffer EOS with K = 263 MeV based on the 
Hartree-Fock approximation, Marek et al. found no explosions 
for the same progenitor model, whereas they indeed obtained an 
explosion for Shen E OS that is even stiffcr with K = 281 MeV 
(see, e.g.. Uankall20T2T ). 



after the shock-stall, on the neutrino-driven shock re- 
vival and the shock propagation for more than 500 ms 
after core bounce. 

The paper opens with the description of the the nu- 
merical method and the employed EOS in Section 2. The 
results of ID and 2D simulations are presented in Sec- 
tion 3 and 4. We summarize our results and discuss their 
implications in Section 5. 

2. SUPERNOVA MODEL 
2.1. Hydrodynamics 
The basic evolution equations are written as follows, 



p— = -VP-pV$ 
dt 



(1) 



(2) 



de* 



P) v] = -pv ■ V$ 
A $ = 4nGp, 



(3) 
(4) 

with rest-mass density p, total pressure P (including neu- 
trino contributions), matter velocity v, total energy den- 
sity e* , gravitational potential $ as well as neutrino heat- 
ing/cooling rates Q v . 4; is the Lagrangian time deriva- 
tive and is the Eulerian time derivative. 

Our 2D simulations are performed using a code which 
is based on the s pectral neutrino transpor t scheme 
IDSA, develope d by ILiebendorfer et al l (12009ft. and the 
ZEU S-2D code (IStone fc Normanlll992tlSuwa et al.ll2010L 
l20TlT) FI Details about the IDSA are given in the next 
subsection. 

In this work we use the conservation form for the to- 
tal energy rather than the internal energy equation used 
in the original ZEUS-2D code to improve the accuracy 
of the total energy conservation. The source term of 
the self- grav ity in Eq. (EB is co mputed by the method 
described in IMiiller et alT(|2010l) . The violation of the 
energy conservation remains within 0.03% of its grav- 
itational binding energy (~ 10 53 erg; see Appendix of 



ISuwa et alJ 120111 ) . The simulations are performed on 
a grid of 300 logarithmically spaced radial zones up 
to 5000 km and 128 equidistant angular zones covering 
< 8 < 7r. For neutrino transport, we use 20 logarith- 
mically spaced energy bins reaching from 3 to 300 MeV. 

Our code is written in spherical coordinates, which 
have the following advantage: We can compare ID and 
2D results of the "same" code. As for ID simulations, we 
just use "one" grid for the lateral direction. This is im- 
possible for a Cartesian grid code or a distor ted grid code 
such as VULCAN/2D ([Burrows et al.l [20061) . In fact, re- 
sults obtained in ID and 2D simulations show exactly the 
same features for the preshocked region (see Figure 1 of 
ISuwa et aLll2010ft . Furt hermore, using the detaile d com- 
parison study of, e.g., ILiebendorfer et ail (|2005j ) based 
on spherically symmetric simulation s, it is possible t o 
validate our code (sec Appendix in ISuwa et alJ 120 111) . 

9 Our code was developed for core-collapse simulation and was 
used for the core co llapse of very massive Population III stars 
l|Suwa c t al. 2007a b, 2009), where a l eakage scheme for ne utrinos 
was implemented l|Kotake et al . 2003; Takiwaki et al. 20Q3). 
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Therefore, we can directly identify the multi-dimensional 
effects by using this code both in ID and 2D simulations. 
On the other hand, the spherical coordinates have a dis- 
advantage for the CFL condition because the length of 
the lateral direction at the center becomes much smaller 
than the length of the radial direction, which can be writ- 
ten as ri A0 - 0.02ri(128/7V e ). Here, n, AO, and N e are 
the mesh width of the radial direction at the center, that 
of the latitude, and the mesh number of the lateral direc- 
tion. This implies that the time step of a 2D simulation 
with this lateral resolution must be shortened by a fac- 
tor of 50 compared to a ID simulation. Therefore, in 
order to perform long-term simulations in spherical co- 
ordinates, we should employ some sort of approximation. 
In this study, the innermost core of ^5 km radius is com- 
p uted in spher i cal sym metry, similar to what was done 
in lBuras et all (|2006a| ). 

2.2. Neutrino transfer 

In this subsection, we describe the scheme of the neu- 
trino transfer used i n this study. This s cheme is based 
on the description in ILiebendorfer et al.l <|2009f ) . 

We solve the transport equation for neutrino distribu- 
tion function fi(t, r, \i, E) depending on the time t, radius 
r, and the momentum phase space spanned by the angle 
cosine fi, of the neutrino propagation direction with re- 
spect to the radius, and the neutrino energy E. In this 
calculation, we include two neutrino species (I = v c ,v c ). 
We divide the distribution function of neutrinos into two 
parts as 

fl = fi+fl, (5) 

where / is an isotropic distribution function of trapped 
particles and f s is a distribution function of streaming 
particles, representing neutrinos of a given species and 
energy which find the local zone opaque or transparent, 
respectively. We solve each component separately as 

D(ft)=ji-(ji + Xi)fi-Xi, (6) 
D(f l s ) = -Ul+Xi)fi+X l , (7) 

where DQ is an operator describing particle propagation, 
j is a particle emissivity, x is a particle absorptivity, and 
£ is the diffusion source, which converts trapped par- 
ticles into streaming particles and vice versa. In this 
app roximation, S is determi ned by the diffusion limit 
(see ILiebendorfer et al~l 120091 for detail) . Weak interac- 
tions are implemented as described in iBruennl (|1985h . 
We assume that the spectrum of /' is thermal and treat 
trapped particles as a fluid element. We characterize the 
thermal equilibrium by a particle number fraction, Y , 
and a particle mean specific energy, Z l , 

where nib is the baryon mass, p is the baryon density, 
h is the Planck constant, c is the speed of light, and 
e u is the energy of neutrinos. The advection of these 
values is calculated in the same way as hydrodynamic 
quantities. Using the net interaction rates, s/ = ji — 
{ji +Xi)(fi between matter and neutrinos, we can 



write the following changes of the electron fraction, Y c 
and the internal energy, e: 

dY c m b Aire 2 f 2 
de mi, Aire 2 f o 

The changes of the electron fraction and specific inter- 
nal energy feed back the emissivity and absorptivity into 
Eq. ©, so that we should find the consistent solution 
by iterating Eqs.(H>|), (fTU)) , and (|TT|) . Once the consis- 
tent solution has been found, Eq. © gives the updated 
distribution function and specific energy of the trapped 
particles, which leads to updated Yi and Zi values by 
Eqs. © and ©. 

2.3. Equation of State 

The equation of state (EOS) in supernova simulations 
needs to handle a large variety of thermodynamic condi- 
tions that reflect different nuclear regimes. For tempera- 
tures below 0.4 MeVF 7 !. time-dependent nuclear reactions 
are important and the baryon composition is dominated 
by heavy nuclei. Above this temperature nuclei are in 
chemical and thermal equilibrium, known as nuclear sta- 
tistical equilibrium (NSE) , where the baryon composition 
is given by temperature, baryon density and proton-to- 
baryon ratio (or equivalently the electron fraction in the 
absence of, e.g., muons). At high densities, typically 
above normal nuclear matter density, the transition to 
a state of matter composed of deconfincd quarks may be 
possible. In addition to the baryons (or quarks), a gas of 
degenerate and possibly relativistic electrons is present. 
Also (non-degenerate) positrons and photons must be 
taken into account. 

In this study we employ in supernova simula- 
tions the two standard EOS for matter i n NSE, the 
ILattimer fe Swestvi (|1991D EOS (LS) and the lShen et all 
dl998j) EOS (SHEN). LS is based on the incompressible 
liquid drop model for nuclei, combined with dripped nu- 
cleoids. The transition to uniform nuclear matter is calcu- 
lated via a Maxwell phase transition. Contributions from 
electrons/positrons and photons are taken into account 
intrinsically. LS has a symmetry energy of 29.3 McV. 
Three different values of the incompressibility parameter 
are available, K = 180 MeV (LS180), 220 MeV (LS200), 
and 375 MeV (LS375). With these one obtains maxi- 
mum neutron star masses of ~ 1 .8Mm (LS180), ~ 2.0M (^ 
(LS220), - 2.7M Q (LS375) fsee lO'Connor fc OttllMTh . 

SHEN gives only contributions for baryons valid for 
matter in NSE. It is based on the relativistic mean field 
theory and the Thomas-Fermi approximation for the de- 
scription of heavy nuclei. The symmetry energy and the 
incompressibility are 36.9 MeV and K = 281 MeV, re- 
spectively. It r esults in a maximum neutron star mass of 
~ 2.2 M m (see lcTCo nnor fc Ottll20TTI : IKiuchi fc KoTakel 
HH). 

Although LS180 fail s to fulfill the re c ent n eutron-star 
mass constraint from iDemorest et al.l (|2010f ) of 1.97 ± 

10 Although the threshold temperature is slightly higher (around 
0.5-0.8 MeV, depending on density and isospin asymmetry), we 
have tested also the higher flash temperatures and found no signif- 
icant difference. 
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Figure 1. Evolution of the central density, comparing the differ- 
ent EOS under investigation (LS180: red solid line, LS220: blue 
dashed line, LS375: green dotted line, and SHEN: gray dot-dashed 
line). 

0.05A/q, we nevertheless include it in our EOS compari- 
son study here because it is widely used and because of its 
similarity to LS220, which does fulfill the Dcmorest et al. 
constraint. 

There are only few experimental constraints on high- 
density and finite-temperature nuclear EOS. For zero- 
temperature symmetric nuclear matter, composed of 
equal numbers of protons and neutrons, the symmetry 
energy can be well determined. Moreover, recent con- 
straints from nuclear theory, experiments and observa- 
tions predict a symmetry energy between 31-33.6 MeV 
with a density gradien t between 49.1 and 80 MeV (see 
lLattimer &; Limll2012t ). These parameters produce cold 
neutron star radii between 11-12 km, in agreement with 
the results of recent efforts to consistently include three- 
body forces ba sed on low-energy chiral perturbation the - 
ory (see, e.g., ISteiner et all 120101 IHebeler et all l2010h . 
Heavy-ion collision experiments can only hardly be used 
to constrain the supernova EOS, because the conditions 
in the collisons are rather isospin-symmetric such that 
they only cover slightly neutron-rich conditions where 
Y e ~ 0.485. But matter in supernova interiors can be 
extremely neutron rich with Y e = 0.01-0.3. Further- 
more, the incompressibility constraint of 200 MeV < 
K < 300 MeV is still weak. 

In this paper, we investigate the impact of the EOS on 
the dynamics of core-collapse supernovae in ID and 2D 
simulations. Since the employed EOS cover a wide range 
of nuclear parameters, we can identify the systematic 
features that are related to the different EOS. 

3. SIMULATIONS IN SPHERICAL SYMMETRY 

In the following subsections we compare results ob- 
tained in simulations using the four EOS which were 
discussed in ^2.31 We u se a 15 M progenitor from 
iWooslev fc Weaver] (|1995|) . We will start with spheri- 
cally symmetric simulations focusing on the differences 
between the different LS EOS. Further below, we will 
also discuss results obtained using SHEN. 

The supernova evolution can be separated into 
the following phases; core collapse, bounce, prompt 
shock propagation, neutrino burst, and accretion phase 
(|Kotake et all 120061: Uanka et al J 120071 ). Note that for 
the post-bounce times considered, exp losions could not 
be obtained in spherical symmetry (see lSuwa "et~aT1[2010l 



[20Tll for more details). 

3.1. Comparison between different LS EOS 

Figure [1] shows the time evolutions of the central den- 
sity. Because of the difference of the incompressibility 
K (large K implies a stiff EOS for otherwise identical 
nuclear parameters), LS180 (red solid line) results in the 
highest central density, while LS375 (green dotted line) 
corresponds to the lowest. However, the central density 
differences between these EOS is < 10%. Since LS220 ex- 
hibits almost identical features as LS180 during the post- 
bounce phase up to about 500 ms and because LSI 80 
has been used in many previous publications of super- 
nova simulations, we mainly compare LS180 and LS375 
in the following. Figure [2] shows the radial bounce pro- 
files of velocity, Y e , temperature and entropy per baryon 
as a function of the enclosed mass. The slightly differ- 
ent temperatures obtained right behind the shock and 
at the very center for the different LS EOS, are due to 
the different incompressibilities. Highest (lowest) tem- 
peratures behind the shock, as well as lowest (highest) 
temperatures at the stellar core, arc reached using LS375 
(LS180) (see Figure [2]). Moreover, the lowest (highest) 
Y e is obtained using LS375 (LS180). However, the dif- 
ferences obtained are small and all LS EOS have very 
similar profiles. Thus, despite the difference of the cen- 
tral density, the structures of the inner cores at bounce 
are very similar for all LS EOS, especially the mass en- 
closed inside the bounce shock and hence the energetics 
of the bounce shock is practically the same. 

All LS EOS show a similar shock evolution post bounce 
up to about 200 ms post bounce, expanding to ~ 250 km. 
After that, the shocks stall and turn into standing accre- 
tion shocks. After about ~ 200 ms post bounce, LS375 
has a larger shock radius than LS180 and LS220 (see 
Figure 3(a) ). 



3.2. Comparison between LS and SHEN EOS 

SHEN exhibits a significantly lower central density 
compared to LS375 which is at first view counter- 
intuitive because SHEN has a lower incompressibility 
(K = 281 MeV) than LS375 (K = 375 MeV). This as- 
pect can be related to the larger adiabatic index of SHEN 
at high densities around w 2 x 10 14 g cm~ 3 . It means 
that there are regimes where curves of constant entropy 
per baryon in a p — p diagram have a steeper gradient 
for SHEN. Figure [4] shows the pressure as a function of 
the density for the employed EOS, at fixed entropy per 
baryon of 1 k^ and electron fraction of Y e = 0.3. Only 
at very high densities, > 3 x 10 14 g cm~ 3 , LS375 has a 
steeper density gradient and becomes stiffer than SHEN, 
mainly due to the higher K. This example illustrates 
that neither incompressibility nor symmetry energy alone 
determine the nuclear properties of the EOS. It is rather 
the entire ensemble of nuclear parameters, and in par- 
ticular also the gradients of the incompressibility and 
symmetry energy, which determine the EOS and hence 
the possible outcome of dynamical simulations. 

At core bounce and shortly after, the shock wave prop- 
agates at high densities, where the highest pressure of 
SHEN results in the largest shock radius. Later during 
the post-bounce evolution, the bounce shock propagates 

g cm" 



towards lower densities on the order of 10 9 



where all EOS are very similar (Figure 3(a) I 
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Figure 2. Radial bounce profiles as a function of the enclosed baryon mass. The velocity jump (at ~ O.7A/0 for LS EOS and ~ O.8A/0 
for SHEN) mark the edge of the inner core, above which the accretion flow is supersonic. The inner core at this time is not swept by the 
shock wave so that the inner core maintains its low entropy. 
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(a) Shock evolution (b) Electron fraction evolution 



Figure 3. The left panel shows the evolution of the shock radius. Note that the time axis uses two different scales to clearly show the 
early post-bounce phase. The right panel shows the evolution of the central electron fraction for the different EOS under investigation. 
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Figure 4. The pressure as function of the density for the case of 
s=lkB/baryon and Y e = 0.3. These conditions correspond to the 
most central part of the collapsing core where the neutrinos are 
completely trapped. Above nuclear density (sa 3 X 10 14 ) g cm -3 
the difference between the EOS is significant. 



Bounce conditions of selected quantities are illustrated 
in Figure [2] in comparison to the LS EOS. The most im- 
portant difference is the larger mass enclosed inside the 
core at bounce for SHEN. This difference cannot only 
be related to the incomprcssibility. More relevant is the 
symmetry energy, which mainly affects the evolution of 
the electron fraction Y e . Note that the dominant con- 
tributions to the pressure are given by the degenerate 
electron gas during the core-collapse phase, while after 
bounce the dominating contributions come from nucle- 
ons in the hi ghly dissociated regime behind the shock 
wave. Figure 3(b) shows the evolution of the central Y e 
for the different EOS. Y e is very similar for all LS EOS. 
But for SHEN, Y e is significantly larger than for any LS 
EOS, already during core collapse. At bounce this dif- 
ference results in Y e ~ 0.32 for SHEN and Y e ~ 0.3 for 
LS. It indicates a slower iron-core dclcptonization dur- 
ing collapse and explains the larger core mass at bounce 
of about 0.8 M Q for SHEN. The models based on LS 
EOS result in core masses of about 0.7 Mq. The dif- 
ferent shock positions at bounce relate to these differ- 
ences in the structures of the core (see the temperature 
and Y e profiles in Figure [5]) and hence reflect the differ- 
ent shock energetics during the initial propagation. Note 
the lower central entropy per baryon for SHEN in com- 
parison to LS, which is mainly due to the well known 
problem of the interactive LS EOS implementation, pro- 
ducing a systematically too low fraction of a-particlcs 
at low density. We did not correct for this in the cur- 
rent study. Note also the lower entropy per baryon for 
SHEN outside about 1.5 Mq. It is related to the different 
non-NSE treatments, based on the ideal Si-gas approx- 
imation, in LS and SHEN. Our findings are consistent 
with previous st u dies c ompa ring LS and SHEN, see e.g. 
iSumivoshi et al.l (|2005t) and IHempel et al.l (|2012[ ). It is 
also worth mentioning that the incorrect alpha-particle 
fraction of LS EOS has a negligible impact on neutrino 
heating/cooling in the post-shock regions because all nu- 
clei (including alpha-particles) diss ociate into free nucle - 
ons (e.g., panel (e) of Figure 7 in IHempel et all 120121) . 
Moreover in the preshock regions, inelastic neutrino scat- 
tering of lig ht nuclei, whic h w as proposed to prehe at ma- 
terial there (jHa xton 198jl see lOhnishi et al.1l2007l for col- 



lective references therein), has recently been shown to 
have neglig ible impact on assisti ng the neutrino-driven 
explosions (lLanganke et al.ll2008ft . 

Note that the inner-core mass of 0.7 and 0.8 Mq 
for the enclosed baryon mass at bounce are too high 
values that are due to the approximations employed 
in our simulations. We omit neutrino scattering on 
electrons/positrons and general- relativistic effects, which 
lead to smaller inner core masses at boun ce (see 
ILiebendorfer et al.|[200lT: iThompson e Fall [2001 ) In ad- 
dition , we use the electron-capture rates from IBruennl 
(|1985f ). also for heavy nuclei, which assume a single 
representative nucleus with average charge and mass. 
Improved electron-cap ture rates on heavy nuclei from 
lLanganke et al.l (|2003l) , which are based on the NSE dis- 
tribution of he avy nuclei, also l ead t o smaller core mas ses 
at bounce (see iHix et all 120031 and lJanka et al.ll2007l for 
details). 

These size differences of the inner core mass lead to 
different shock evolutions. Due to the larger inner core, 
SHEN initially shows a slightly more rapid expansion of 
the shock wave. However, ~ 50 ms after bounce the 
shock contracts faster towards the PNS for SHEN than 
for any LS EOS. The reason is a larger energy loss with 
SHEN when the shock propagates with its larger energy 
across the neutrinosphcrcs, from where the deleptoniza- 
tion burst is released. This effect can also be seen in 
the luminosities. Before the explosion is launched, the 
electron (anti)neutrino luminosity is determined by mass 
accretion during the post-bounce evolution. With the 
accretion rate M and an enclosed mass M at the neutri- 
nospheres R v one obtains the accretion luminosity as 



L acc ~8xl0 52 



M 



1.5 Mi 



50 W 



M 



t erg s 

1 Mrs" 1 ) 5 



Q 2 ) 

Figure [5] shows the post-bounce evolution of this quantity 
for the different EOS under investigation. The highest 
mass accretion rates for SHEN explain the highest lu- 
minosities for SHEN until about 150 ms post bounce, 
in comparison to the LS EOS. During the later post- 
bounce evolution, the mass accretion rates at the neutri- 
nospheres become very similar for all EOS considered. 
Only small differences can be identified with slightly 
lower luminosities for the stiffcr LS375 and SHEN on the 
one hand, and slightly higher luminosities for the softer 
LS180 and LS220 on the other hand. 

4. RESULTS OF 2D SIMULATIONS 

In this section, we present our results of 2D simula- 
tions. We perform 2D simulations for LS180, LS375, and 
SHEN. Since LS220 shows almost the same features as 
LS180 as seen in the previous section, we do not include 
it in this section. 

4.1. Comparison between ID and 2D using LS180 

First of all, we compare the results of spheri- 
cally symmetric and axially symmetric simulations us- 
ing LS180. For the 2D simulation we obtain a 
weak explosion , corresponding to the so-called "p as- 
sive expansion" ()Buras et al. 2006a: S uwa et aLll2010f) , in 
which the rate of mass accretion through the shock is 
larger than that of mass-ejection, so that the mass of the 
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Figure 5. Post-bounce evolution of the accretion luminosity for 
the EOS under investigation. 
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Figure 6. The mass trajectory as a function of time after the 
bounce for 2D simulation of 15 Mq with LS180. The thick black 
line represents the angle-averaged shock radius. The grey lines 
show the mass from 1.0 to 1.7 Mq at intervals of 0.01 Mq. Two 
thin black lines indicate 1.4 and 1.5 Mq, respectively. 

central PNS is still increasing after the onset of shock 
expansion^ 1 " 1 !. On average, only the continuous shock ex- 
pansion to larger radii is obtained with yet no significant 
matter outflow (see Figure [5]) for the simulation times 
considered. The shock expansion is an effect of the accu- 
mulation of hot matter in combination with a decreasing 
ram pressure of the material ahead of the shock. The ram 
pressure decreases because of a decreasing accretion rate 
and the propagation of the shock to larger radii where 
the gravitational potential is less deep. 

In Figure [7] we show the entropy distribution at the 
north pole (top panel) and the south pole (bottom panel) 
for the 2D simulation. The deviation from a spherically 
symmetric distribution becomes prominent after about 
20 ms post bounce when convection starts behind the 
shock and shock wave oscillations start to introduce sub- 
stantial differences of the shock position at the north 
and south poles. The standing accretion shock insta- 
bility (SASI) drives such shock oscillations on a longer 
timescale on the order of several 100 ms. Until about 
300 ms post bounce, the standing accretion shock ex- 
pands only slowly to several 100 km. After that, the 

11 Note that the post-shock material have partly positive veloc- 
ities (e.g., Figure 9(a)). 



shock expands faster and reaches about 1000 km at the 
north pole and 2000 km at the south pole at about 500 ms 
post bounce. We will return to the origin of the increas- 
ing difference between north- and south-pole evolution 
further below. Here, we will discuss the entropy distri- 
bution obtained in comparison with the spherically sym- 
metric simulation. 

Figure [8J shows the radial entropy distribution ob- 
tained in the 2D simulation (red and blue points) and 
the ID simulation (green line), at 50 ms, 100 ms, 
200 ms, and 300 ms after bounce. In general, the en- 
tropy spreads over a wider range in 2D than in ID, 
because of the convective activity between the shock 
and the gain radius (i.e., the region where neutrino 
heating exceeds cooling). The higher entropy in 2D 
is realized by more efficient neutrino heating, in com- 
parison to ID (iHerant et all 119941 IBurrows et al.l [19951: 
iJanka fc Muellerlll996l; IMurphv k Burrowsll2008l) . Note 
that at late times, shown in the bottom panels of Fig- 
ure [8j once maximum entropy (s ~ 13 ks per baryon in 
the bottom left panel and 16 in the bottom right panel) 
is obtained, matter at small radii (~ 100 km) does not 
gain additional entropy/encrgy as it floats to larger radii. 
This is related to the convection which drives matter 
from the gain radius, where the heating is maximum, to 
larger radii, where the heating is less efficient. Thus, the 
specific entropy is limited to the value at the gain radius. 
This is illustrated in Figure [FJ where red points represent 
matter with positive radial velocity (i.e., escaping) that 
have high entropy, while blue points (accreting material) 
have low entropy. For comparison, we show the weak 
equilibrium entropy for the ID simulation (green dot- 
ted line) , that would be achieved by infinitely- long expo- 
sure of matter to the prevailing neutrino abundances and 
spectra. We construct this line by keeping hydrodynamic 
quantities unchanged and solve the neutrino emission 
and captures until equilibriu m is almost achieved (see 
also iLiebendorfer et al.ll2004l) . The region where the en- 
tropy is below this line has the potential to be heated by 
neutrino radiation up to the equilibrium entropy, while 
the region above this line is dominated by neutrino cool- 
ing. In ID the entropy beyond the gain radius (corre- 
sponding to the radius of the maximum entropy, r ~ 100 
km) decreases with increasing radius, while the equilib- 
rium entropy increases continuously with radius. The dif- 
ference between the obtained entropy in the simulations 
and the equilibrium entropy is caused by the insufficient 
time for the heating, i.e., too short advection timescale 
compared to the neutrino heat ing timescale (see the dis- 
cussion in iBuras et all l2006al for instance). Thus, the 
region above the gain radius has room to absorb more 
neutrinos in multi-dimensional simulations. The convec- 
tive motion continuously exchanges equilibrated matter 
from the gain radius by below-equilibrium matter from 
layers that are too distant for efficient neutrino heat- 
ing. This exchange of matter by fluid instabilities en- 
ables continued neutrino heating at the gain radius in 
combination with the stead y distribution of th e heated 
matter to the outer layers (Herant et al.|[l994T) . More- 
over, the convective motion is expected to prolong the 
advection timescale and lead to more effective neutrino 
heating. In summary, the higher entropy obtained in 2D 
simulation is due to the presence of convection and SASI 
in the gain region. However, even in 2D the advection 
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Figure 7. Post-bounce evolution of the entropy distribution, illustrated in a space-time diagram. The top and bottom panels correspond 
to the northern pole and southern pole, respectively. The high-entropy region (red or yellow) shows the post shock layer. The shock wave 
oscillation produced by SASI and convective motion is apparent. The shock wave at the southern pole reaches a distance of 2000 km from 
the center earlier than at the north pole, corresponding to the existence of a I = 1 mode of the SASI. 



timcscalc is finite so that the equilibrium entropy is not 
realized. Part of the material just below the gain radius 
maintains a higher entropy than the equilibrium value 
because of the overshooting of high entropy matter by 
convective motion. This phenomenon disappears deeper 
inside the gain region (r < 50 km) because the weak in- 
teraction rates at high density become large so that the 
entropy coincides with the value determined by weak and 
thermal equ ilibriu m. 



In Figure 9(a) we show the time snapshot of the en- 
tropy distribution at 400 ms after bounce. The left panel 
shows the central part and the right panel is a wider 
view. In the left panel, red arrows indicate infalling mat- 
ter (negative radial velocity) and blue arrows outstream- 
ing matter (positive radial velocity). An important 
structure becomes visible: low-entropy infalling matter 



( "downflow" , see 


Herant et al.lll992HBurrows et al.lll995t 


Uanka & Mueller! 


19961) is accreting onto the central PNS 



as indicated by the green area in Figure 9(a) (left panel). 
On the right-hand-side panel of Figure |9(a)[ a triple 
shock point is seen to be located in the northern hemi- 
sphere near the position of (600 km, 200km). If one now 
compares with the left- hand-side panel of Fig 9a, one rec- 
ognizes that the cold downflow through the triple point 
advects directly to the surface of the PNS. Hence, the 
cold downflow stems from large distances, between 500 
and 1000 km, and reaches down to the PNS surface. It 
exists as a quasi-stationary structure and supplies en- 
ergy and momentum to the PNS in an asymmetric way. 
The momentum is reflected at the surface of the PNS 
and becomes visible as a pressure wave in Figure 1131 
which pushes the shock wave to larger radii and increases 
the amplitude o f the SASI by the s o-called "advective- 
acoustic cycle" (jFoglizzo et al.l 120071 ). This feedback cy- 
cle stands at the origin of the anisotropic shock evolution 
between north and south poles which significantly pro- 
motes the shock expansion reaching larger radii at the 
south pole. 

The IDSA neutrino transport scheme employed in this 
study permits an estimate of the proton-to-baryon ra- 
tio of the possibly ejected material from the PNS sur- 
face. In Figure ITOl we show the Y e distribution at 400 ms 



post bounce. Infalling material at large radii, domi- 
nated by heavy nuclei, starts with the progenitor value 
of Yg ~ 0.5. The electron fraction then decreases contin- 
uously due to electron captures as the material is com- 
pressed while accreting onto the PNS surface, where Y e 
finally reaches values as low as ~ 0.05-0.2. However, 
the high-entropy region, which corresponds to expand- 
ing matter, experiences a continuous flux of neutrinos 
that stream off the PNS surface. There, in the presence 
of similar electron neutrino and antineutrino luminosi- 
ties, and small average energy differences below 4 MeV 
between u e and D e , matter becomes proton rich. Similar 
results have been reported in previous studies of massive 
star explos ions that employ sophisticated neutrino trans 
port (e.s... iLiebendqrfer et al.l 120011 iBuras et al.l 12006b 



Frohlich et al.1 120061 : iPruet et all 120061: iWanaiol 12006 : 
Fischer et al.ll2010D . 

In contrast to spherically symmetric models, multi- 
dimensional models permit mass outflow alongside ac- 
creting matter and hence, e.g., mixing with the sur- 
rounding material. Mixing opens interesting aspects for 
the chemical evolution of explosion models which cannot 
be obtained using ID models. This has been discussed 
at the e xample of the 8 .8 Mp O-Ne-Mg-core explosion 
model in Wa naio et al.l (|2011[ ) by comparing ID and 2D 
results. Neutron-rich conditions were obtained in both 
cases due to the extremely fast initial expansion of the 
ejecta. However, while only slightly neutron-rich condi- 
tions could be obtained in ID, mushroom-like and more 
neutron-rich pockets developed in 2D. 

4.2. Comparison of 2D models with different EOS 

Here, we compare the 2D results for the 15 Mq pro- 
genitor using the different EOS, i.e., LS180, LS375, and 
SHEN. Figure [11] shows the shock-radius evolution for 
each model. Each of which has three lines, corresponding 
to maximum, angular averaged and minimum shock ra- 
dius (from top to bottom) . Note that the minimum shock 
radius often coincides with the position of the triple point 
of the shock, whi ch is t he starting point of the cold down- 
flow (see Figure 9(a) I. In Figure HT1 we can see that 
LS180 and LS375 (red solid lines and green dotted lines) 
show similar trajectories. There is a gradual expansion 
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Figure 8. Radial entropy distribution of the simulations with LS180. Red and blue points represent result of the 2D simulation, the green 
solid line corresponds to the ID results. In addition, the green dotted line shows the entropy in weak equilibrium that would be achieved 
by infinitely-long exposure of the matter to the prevailing neutrino abundances and spectra. 



phase during the post-bounce time i p b < 200-300 ms, 
where the maximum shock radii grow to about 500 km. 
While afterwards, they grow more rapidly in both mod- 
els and reach about 2000 km for LS180 and 1800 km for 
LS375 at 500 ms post bounce. These two models show 
no phase of shock retraction. For the 2D simulation us- 
ing SHEN, on the other hand, the shock wave does not 
enter this second phase of fast shock expansion even at 
late times (see Figure [TTj) . 

Before analyzing the the origin of this observed dif- 
ference in the shock evolution between SHEN and LS180 
further, a comment regarding the likelihood to find explo- 
sions in our models seems in order: At the time the 2D 
simulations are sto pped, the "diagn osing explosion en- 
ergy", as defined in lSuwa et al.ll2010l achieves ~ 10 50 erg 
in both LS EOS cases, while explosions for the 2D simula- 
tions using SHEN are unlikely to occur. LS180 has been 
applied in axially symmetric simu l ations of the same 15 
Mq progenitor in lMarek fc Jankal ()2009l ). where the av- 
erage shock radius was found to oscillate around 200 km 
for several 100 ms until it eventually starts to expand at 
about 600 ms post bounce, i.e. much later than in our 
simulations. The fast and continuous shock expansion 
to increasingly larger radii for our LS180 simulation is 



likely due to the omission of the continuous energy loss 
by the emission of heavy lepton neutrinos, and possibly 
the neglect of general relativistic effects in the gravita- 
tional potential. Coming back to the discussion of the 
differences between LS180 and SHEN, we will analyze 
in the following the two different aspects of (a) neutrino 
heating/cooling, and (b) subsequent convective activity. 

Figure [12] illustrates the evolution of the neutrino lu- 
minosities and average energies in ID (thin lines) and 
2D (thick lines), comparing LS180 (red solid line) and 
SHEN (grey dot-dashed line) , for electron neutrinos (left 
panel) and antineutrinos (right panel) , respectively. Note 
the rapid variations in the luminosities on a millisecond 
timescale, which are due to the convective a ctivity in the 
vicinity of the PNS (e.g. JMarek et al.ll2009D . At first, the 
luminosities (sampled at the equator) are in qualitative 
agreement with the results obtained in spherical symme- 
try. The mean neutrino energies in the 2D models are 
slightly, but systematically, smaller than in the ID mod- 
els. When the standing shock starts to propagate faster 
to increasingly larger radii in 2D, i.e. at about 350 ms 
post bounce for LS180, the v e and v e energ ies cease t o 
rise and stay nearly constant (see Figs. 12(a) and 12(b) I. 
Furthermore, the luminosities decrease rapidly and dc- 
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Figure 9. Snapshot at 400 ms post bounce of the entropy distributions comparing LS180 and SHEN. Red and blue arrows show the in- 
and outflowing velocity fields, respectively. The left panel shows the inner part while the right panel is a wider view. The shock wave is 
strongly deformed due to SASI activity and turbulent matter motion induced by convection. In the left panel we can see a cold downflow 
of matter that comes from the northern hemisphere, which activates a pressure wave in the southern hemisphere (see text for detail). This 
feature originates from the triple point of the shock wave that is shown in the northern hemisphere of the right panel. 
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Figure 10. Y e distribution for the simulation using LS180 at 
400 ms post bounce. Escaping material exhibits proton-rich condi- 
tions with Y e > 0.5, while infalling material becomes neutron rich 
due to electron captures. 



viate from the ID results (see Figs. 12(c) and 12(d)). 
This phenomenon is more prominent tor v e , whose neu- 
trinosphcrc is located at lower densities than those of u e . 
The LS180 and SHEN EOS lead to different central 



Figure 11. Shock trajectories for the 2D simulations using the 
different EOS LS180 (red-solid lines), LS375 (green-dotted lines), 
and SHEN (gray dot-dashed lines). We show three lines for each 
model, which correspond to the maximum shock radius, the angu- 
lar averaged one, and the minimum radius, from top to bottom. 

configurations during the post-bounce evolution, as a 
result of the different PNS contraction behaviors. In 
comparison to the stiffer SHEN EOS, the soft LS180 
EOS results in a more compact PNS with higher cen- 
tral density and a steeper density gradient at its surface. 
The impact of different PNS structures on the proper- 
ties of emitted neutrinos has bee n discussed in man y 
places of the literature (see, e.g., iHempel et al.|[201^) . 
A direct consequence of the different PNS structures is 
the d ifferen t ev olution of the neutrino luminosities (see 
Figs. 12(c) and |12(d)[ ). In ID, the neutrino luminosi- 
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ties (thin lines) become slightly smaller for SHEN com- 
pared to LS180. However, the difference between SHEN 
and LS180 is much sma l ler than what has be en found in 
iSumivoshi et~aTl (pOOl : iFischer et afl (pOOl . The rea- 
son might be that our PNS is less compact due to the 
neglect of heavy-lepton neutrinos and general rclativistic 
gravity. As the shock for LS180 in the 2D simulations 
starts to expand continuously after about 350 ms post 
bounce, we will focus on post-bounce times before that 
moment to search for the origin of more efficient neutrino 
heati ng in LS1 80 com pared to SHEN. It can be seen in 
Figs. 12(c) and 12(d) that the v e luminosities are rather 
similar lor LSI 80 and SHEN, while the v e luminosities 
differ considerably during the first 250 ms post bounce. 
We explain these differences by the fact that the v e neu- 
trinosphcres are located at higher density where the EOS 
differences between LSI 80 and SHEN are larger than at 
the lower density where the v e neutrinospheres are lo- 
cated. The sizeable difference in the v e luminosity leads 
to less efficient heating for SHEN, and consequently to 
a lower entropy per baryon in the gain re gion t han in 
the simulations based for LS180 (see Figure [9 (b)| . How 
the luminosity decays at later times, as observed here 
after the onset of the explosion in the simulations using 
LS180, and as predicted by recent spherically symmetric 
explosion models, needs to be elaborated by long-term 
multi-dimensional simulations that cover several seconds 
after the onset of explosion. This is yet computation- 
ally very expensive and beyond the scope of the present 
study. 

In the following paragraph, we will return to aspect 
(b), which is a comparison of convective activities be- 
tween our 2D models using SHEN and LS EOS. The left 
panel of Figure [13] shows the post-bounce evolution of the 
kinetic energy of the lateral (9) direction, J pvg/2 dV, 
where V is the entire simulation domain. This quantity 
is zero in spherical symmetry such that non-vanishing 
values imply a deviation from sphericity. Models based 
on LS180 and LS375 exhibit a significantly larger am- 
plitude than SHEN for i p b < 50 ms. This is the so- 
called prompt convection, which is triggered by the neg- 
ative entropy gradient left behind the prompt shock as it 
weakens by dissociation of nuclei and eventually comes 
to a complete stall by the emission of the deleptoniza- 
tion burst. In the right panel of Figure [T3] we show the 
radial profiles of the root- mean-square of the radial veloc- 
ity, (/ v^svuOdO) 1 ! 2 . For LS180 (upper part), convective 
motion is apparent between 30-100 km during the first 
50 ms after bounce, which is seen as a yellowish island in 
the plot. This shows that the convective flows produced 
by the fluid instabilities advect down to the inner surface 
of the PNS (white dashed line). The occurrence of this 
advection- me diated convective m o tion has been alread y 
mentioned in lBuras et a l. ( 2006a!) ; [Dessart et al.l (|2006| ). 
For SHEN (lower panel of Figure fTBT b)). such a behavior 
is only weakly visible due to the smaller volume of the 
convectively unstable region. Note that after the convec- 
tion has ceased, both, LS180 and LS375, still have larger 
kinetic energy than SHEN. 

In order to analyze the evolution of the shock wave, we 
conduct a mode analysis as follows: The deformation of 
the shock surface is decomposed into spherical harmonic 



components 
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Note that due to axisymmetry, only m = harmonics 
are relevant. The coefficients, ag, can be calculated by 
the orthogonality of the Legcndrc polynomials, 
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The position of the shock surface, R s (0), is estimated 
from the iscntropic surface of ,s = 6fce per baryon. 

Figure [14] shows the post-bounce evolution of the low 
mode {t = 1 and 2) for LS180 and SHEN. During the 
early phase (i p b <100 ms after bounce), LS180 has a sig- 
nificantly larger amplitude than SHEN. After that the 
amplitudes are similar in the time window 100 < i p b < 
300 ms. For LS180, the shock wave begins to expand 
more rapidly at i p b ~ 300 ms, after which the ampli- 
tude is decreasing continuously with time. For SHEN, 
the shock wave continues to oscillate at smaller radii and 
does not expand. This in turn results in an increasing 
amplitude. At < p b > 500 ms, LS180 shows a dominant 
t = 1 mode (red line) due to the expansion of the shock, 
while SHEN shows a decreasing amplitude because the 
retracting shock wave is approaching a spherical shape. 
This corroborates the failed explosion for SHEN. Note 
that the mode analysis of the shock radius, albeit very 
helpful to extract the information of SASI modes, serves 
as a reliable indicator of the explosion only if it is com- 
bined with other criteria such as the convention timescale 
argument regarding the heating/advection timescale and 
the antesonic condition that we will discuss later in this 
section. 

Next, we discuss the EOS dependence of the SASI ac- 
tivity (acoustic- vorticity cycle). We show the pressure 
deviation for LS180 and SHEN in Figure IT5l It is deter- 
mined as follows: 



{Uo [M{r,6)-M(r)] 2 sin effl} 
M{r) 



1/2 



where 



M{r,9) = p(r,6)v*(r,6)+P(r, 9), 
2 



M{r) 



M(r, 9) sin 



(15) 

(16) 
(17) 



M. is the total pressure, including the ram pressure of the 
infalling material. It should assume the same value ahead 
and behind the shock wave (corresponding to the mo- 
mentum part of the Rankinc-Hugoniot equation). There- 
fore, M. is a useful quantity to investigate the shock ex- 
pansion: If Ai is greater behind the shock than ahead 
of the shock, the shock wave propagates outward, and 
vice versa. Before 100 ms postbounce, the dispersion 
of M. is small (~ 0.1) (especially for SHEN; see Fig- 
ure ll5[) . Later, the pressure perturbation increases in the 
shocked region and the shock wave is pushed outward. 
For LS180, the pressure perturbation continues to grow 
after t p b ~ 200-300 ms and a strong pressure wave trans- 
fers momentum to the shock. This phenomenon is absent 
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Figure 12. Post-bounce evolution of neutrino mean energies and luminosities, for the 2D simulation using LS180 and SHEN sampled at 
the equator (thick lines). For comparison, we also show results from the ID simulation using the same EOS (thin lines). 

in the simulations using SHEN. As a consequence in the 
LS180 model, the pressure behind the shock wave be- 
comes greater than the ram pressure ahead of the shock, 
so that the shock begins to propagate outward without 
receding. Note that the pressure wave appears to be pro- 
duced close to the PNS surface (see the thick dashed lines 
in Figure ITS"]) , where the vorticity is reflected at the steep 
density gradient (see also the right panel of Figure [T3")) . 

In order to characterize the spatial anisotropy of down- 
flows for the different EOS used, we estimate the follow- 
ing ratio, 1Z, which is determined by 



n 



El 



1. 



(18) 



where E% and E S K are the kinetic energy accreted from 
the northern (TV) and southern (S) hemispheres, respec- 
tively. Ek = J dA, where the integration extends 
over the surface of a sphere with a radius of 100 km. 
In Figure 16(a)| we show the post-bounce evolution of 
11, comparing LSI 80, LS375, and SHEN. It implies that 
LS180 and LS375 have greatly asymmetric mass accre- 
tion after 300-400 ms post bounce, while for SHEN the 
mass-accretion asymmetry is small. 
In the following paragraphs, we investigate possible 



indicators of the onset of explosion which have been 
previously explored in the literature. One of those is 
the mass-weighted average entropy in the gain region, 
discussed as an indicator of the n eutr ino-heating effi- 
ciency m iMurphv fc Burrows! (|2008l ) and lNordhaus et al.1 
(|2010f K It has been suggested that exploding models ex- 
hibit higher entropy in the gain region prior to the on- 
set of explosion, as a result of more efficient neutrino 
heating, in comparison to non-e xploding models of the 
same progenitor. Figure |16(b)| shows the post-bounce 
evolution of the average entropy, (s), comparing LS180, 
LS375, and SHEN. In agreement with previous studies, 
we find that the optimistic simulations, which show a 
continuous shock expansion to increasingly large radii 
(LS180, LS375), have higher average entropy than the 
less-optimistic one (SHEN), until the onset of explosion 
at about 350 ms for LS180. Once the explosion has been 
launched, the entropy-growth rate slows down and even 
decreases as a consequence of a continuous accretion of 
cold matter from the upwind side. Note that the en- 
tropy per baryon will rise again at the moment when 
mass accretion will turn into mass outflow, such that net 
neutrino heating establishes at the surface of the PNS. 
A similar phenomenon has been observed when increas- 
ing the dimensionality in simulations of neutrino-driven 
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(b) Distribution of the root-mean-square of the radial velocity in a 
space-time diagram. The top half and bottom half correspond to LS180 
and SHEN, respectively. 



Figure 13. Post bounce evolution of kinetic energy and root-mean-square of the radial velocity field as indicators for convective activity. 
In graph (a), both, LS180 (red solid line) and LS375 (green dashed line), exhibit a strong activity, while the activity for SHEN (gray 
dash-dotted line) is weaker by a factor of ~ 3 during the first 50 ms after bounce. After that, prompt convection vanishes and LS180 and 
LS375 maintain slightly larger kinetic energy than SHEN. In graph (b), the colored region is the post-shock layer. The black dash-dotted 
lines indicate the position of the gain radius and the white dashed lines are the radii of p = 10 11 (outer) and 10 13 g cm -3 (inner), where p 
is the angle averaged density. The advection-mediated convective motion (see text) is apparent for 10 11 g cm~ 3 < p < 
tpb < 50 ms for LS180. After such convective motions vanish, weak continued PNS activity is visible above p = 10 13 g cm~ 3 . 
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Figure 14. Post-bounce evolution of the amplitude 1=1 and 2 
SASI modes, comparing LS180 (thick lines) and SHEN (thin lines). 
During the early phase (t p b < 100 ms) LS180 has a larger ampli- 
tude due to the strong prompt-convection activity. Later, LS180 
exhibits a smaller amplitude due to the continuously expanding 
shock wave to larger radii, while SHEN has still an increasing 
amplitude. However, also for SHEN, it saturates and begins to 
decrease at t p b > 400 ms, indicating a failed explosion. 

supernova explosions, e.g., from a xial symmetry to three 
dimensions in par ametric studies (jNordhaus et al.1l2010t 
lHanke et al.ll2012h . The presence of convection and SASI 
exposes matter in the gain region to net neutrino heating 
for a longer period. In our models based on LS180, the 
stronger convective activity leads to more efficient neu- 
trino heating for matter in the gain region such that a 
higher entropy per baryon is achieved than in the models 
based on SHEN. 

Another important indicator to predict neutrino- 
driven explosio ns is the ma ss enclosed i n the gain re- 
gion (see, e.g., Uankall200l . In Fi gure 
that it is always smaller for SHEN than 
LS375 



with lower central density and a less steep density gradi- 
ent at the PNS surface than LS. The mass enclosed inside 
the gain region reaches its maximum at about 100 ms 
post bounce, about 0.065 M Q for SHEN and 0.085 M Q 
for LS. This time corresponds to the time with the max- 
imum accretion rates (see the accretion luminosities in 
Figure [5]). After that, the mass enclosed inside the gain 
region decreases. After about 150 ms post bounce the en- 
closed mass in the gain region for SHEN falls significantly 
below the corresponding values for LS. Note that other- 
wise the evolutions for SHEN and LS are rather similar 
up to this moment, see Figure [T"5l but about 100 ms later 
they differ substantially. As discussed above, this is a 
consequence of more efficient neutrino heating for LS due 
to a higher D e luminosity on the one hand and a larger 
mass enclosed inside the gain region on the other hand. 
For SHEN, the enclosed mass inside the gain region drops 
below 0.02 Mq between 200-300 ms after bounce, while 
it t urns a round at a value of 0.04 Mq for LS (see Fig- 



16(c)) 



16(c) 



find 



0FTSI8O and 

because SHEN produces a less compact object 



As the shock starts to expand for LS only, 
the gam-region mass differences between SHEN and LS 
diverge to an extent that it is not meaningful to compare 
them at a later postbounce time. In summary, the differ- 
ent structures of the protoneutorn stars for SHEN and 
LS lead to a smaller mass enclosed inside the gain region 
and lower luminosities (mainly v e ) for SHEN. Both ef- 
fects inhibit the efficiency of the neutrino heating in the 
models based on SHEN, that do not lead to an explosion 
in our 2D models. 

Additionally the less efficient heating for SHEN can 
also be related to the absence of large-scale hydro- 
dynamic instabilities. Deformations of the standing- 
accretion shock assisted by buoyant convections are sig- 
nificantly s maller for SHEN compared to LSI 80 (com- 
pare Figure 9(a) and 9(b)). Moreover, we cannot find 
any large down-streams of material for SHEN, which for 
LS180 supply continuously energy to the central regions 
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Figure 15. Post bounce evolution of the dispersion of M (Eq. 1156 for LS180 (top panel) and SHEN (bottom panel). The white color is 
set below the minimum value of 0.1. The thick dashed lines show the position of the neutrinosphere representing the surfaces of the PNS. 

in an asymmetric way. 

Because the structure of the PNS now appeared several 
times to stand at the origin between runs using LS180 
and SHEN, we wonder whether the radius of the PNS 
might serve as a reliable explosion indicator. The post- 
bo unce ev olution of the PNS radius is shown in Fig- 
After the rapid initial rise up to about 90 km 



16(d) 



for the models using LS180 and LS375, between 10-20 ms 
post bounce, the PNS radii contract rapidly to about 
80 km at about 50 ms post bounce. Using SHEN, the 
PNS has a smaller maximum radius of 85 km at about 
20 ms post bounce. Later, the PNS contracts on a longer 
timescale on the order of several 100 ms for all EOS, 
where initially SHEN leads to the smallest PNS radius 
up to about 200 ms post bounce. After that, the PNS 
using LS180 contracts fastest and leads to the smallest 
radius. LS375 has the largest PNS radius during the en- 
tire post-bounce evolution. Thus, the PNS radius, which 
is often used as an indicator of the energy budget avail- 
able from core collapse, is not a good predictor for the 
explosion. The same is equally true for the compactness 
(GM m /c 2 R m ). 

In stead, we conjecture that the time derivative of PNS 
radius would be more important, i.e., the faster contrac- 
tion of PNS is better for the shock revival . This is con- 
sistent with the discussion of iJankal (|2012f) . in which he 
argued that the smaller PNS radius is better, not faster 
contraction. It should be noted that in his Figure 4 one 
finds that the model with smaller PNS radius is identical 
to the faster contracting model. Therefore, our conjec- 
ture (faster contraction is better for explosion) is valid 
for their simulations. 

In Figure 1 1 71 we show the a n te-son ic condition pro- 
posed by iPejcha fc Thompsonl (|2012t ). who suggested 
that if max(Cg/Wo SC ) > 0.20 where c s is the sound speed 
and v esc is the escape velocity (= \/2$), the stalled shock 
will turn to an actively expanding shock front. For the LS 
EOS we find a rapid increase of max(cg/z^ sc ) after about 
300 ms post bounce, which corresponds well with the on- 
set of rapid shock expansion. For SHEN, max(Cg/i>o SC ) 
oscillates around its maximum value of about < 0.25 
on a longer timescale on the order of several 100 ms. 

12 This critical value depends on the microphysics (O. Pejcha, 
private communication). 



Note that the strong spike of this indicator observed for 
SHEN originates from the multi-dimensional fluid mo- 
tions (e.g., shock collision with the PNS surface) , which is 
not co nsidered in the ID analysis bv lPeicha fc Thompsonl 
(|20 12f ) . The spikes, stemming from a transient local con- 
dition, disappear soon by energy redistribution by neu- 
trinos as well as by convective motions. Thus, our result 
suggests that transient hot spots of the antesonic con- 
dition do not suffice to predict an explosion. In order 
to have predictive power the antesonic condition should 
be satisfied in a large volume as it is the case in our 2D 
models based on LS EOS . 

The ratio of the advection timescale (r a d v ) and heat- 
ing timescale (rh ea t) is also a well-discussed indicator of 
the explosion (e.g.. iThompson et al.l 120051 : iBuras et al.l 
I2006al iMarek fe Jankal 120091 ) . Thisratio is depicted in 
Figure 1181 If this indicator is greater than unity, i.e., 
T'adv/'Hicat ^ 1, the neutrino heating proceeds fast enough 
to gravitationally unbind the fluid element, otherwise the 
matter is swallowed by the cooling region below the gain 
radius before being heated up. Note that r a d v is deter- 
mined by the travelling time of mass shells between shock 
and gain radii. This is an "effective advection timescale" . 
In 2D simulation, the velocity distribution is stochastic 
because of the turbulent motion and the S ASI activity so 
that the simple estimate by the velocity profile and radius 
between shock and gain radii is misleading. Considering 
a certain mass shell with enclosed mass Mi(r) = J pdV r , 
where V r is the volume of a sphere with a radius of r, 
we note the transit time at the shock and gain radii as 
tshock(Mi) = i(-Mj(r s hock)) and t gailL (Mi) = t(Mi(r ga , in )), 
where r s h oc k and r ga i n mark the shock and gain radii at 
some time. Then, we determine the advection timescale 
as r adv (A/j) ee tg ain (M i ) - t shock (Mi). The angular aver- 
aged values for the shock and gain radii are used. More- 
over, Theat is estimated by the mass weighted integral 
of ebind/Qiy in the gain region, where ebind is the local 
specific binding energy (the sum of internal, kinetic and 
gravitational energies) and Q v is the specific heating rate 
by neutrinos. Figure [18] implies that the models with 
continuous and rapid shock expansion (LS180, LS375) 
exhibit a more rapid increase of the ratio. For SHEN, 
the ratio shows also large oscillations with several peaks 
(even as large as 100). However, this model represents 
the least optimistic case for possible explosions and hence 
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(a) Accretion ratio determined by Eq. JTHJ. 
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(b) Mass-weighted average entropy per baryon in the gain region. 
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(c) Mass enclosed in the gain region. 
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Figure 16. Post-bounce evolution of selected quantities. In graph (a), the models using the LS EOS have a larger anisotropic downfiow, 
i.e. the ratio becomes larger than 100, while SHEN lies between 0.1 and 10 so that the downfiow is rather spherical, slightly oscillating 
from north pole to south pole and vice versa on a timescale on the order of fOO ms. 
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Figure 17. Maximum ratio of the sound speed and the squared 
escape velocity as a function of time after bounce. The models 
with a continuously expanding shock (LSf80, LS375) satisfy the 
criterium max(Cg /v^ BC ) > 0.2, while SHEN shows a very low value 
with the exception of few transient spikes. The critical value in 
this case should be raised to < 0.25. 

the criteria of the timescale ratio is not a robust indicator 
for a successful explosion. 
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Finally, the heating efficiency shown in Figure 19(a) 
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Figure 18. Post-bounce evolution of the ratio of advection and 
heating timescales. The optimistic models (LSf80, LS375) imply a 
much faster heating than advection through the gain region, while 
SHEN exhibits several oscillations, which correspond to the shock 
oscillation and the prolonged dwell time in the gain region. 



is used as a measure for the neutrino he ating process 
(|Marek fc Jankal 120091 ; iMiiller et al.ll2012a|) . The dchni- 
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tion of the heating efficiency is given by 

Qhcat 



?7heat 



(19) 



where Qhcat is the volume-integrated neutrino heating 
rate in the gain region. We did not find any signifi- 
cant difference for this indicator between models with 
different EOS. In addition, we show the time evolu- 
tion of the growth parameter of convection, Xi which 
is used as an indicator for the on set of the convection 
induced by the entropy grad i ent (iFoglizzo et al.l [2006t 
iFernandez fe Thompson! [20091 ; IMiiller et al.ll2012b|) . de- 



fined as 



X 



wbv 



dr, 



(20) 



where lobv is the Brunt- Vaisala frequency. 
iFoglizzo et al.l (|2006l ) showed that x must exceed 



3 for convective motion to glow. Figure 19(b) exhibits 
the time evolution of \ for investigated models and 
suggests that EOS docs not significantly affect the onset 
of the convection induced by the neutrino heating. 

At this point, we remind the reader again that in this 
study we omit important cooling and thermalization con- 
tributions, i.e. heavy-lepton neutrinos and scattering on 
electrons/positrons. However, these simplifications affect 
the evolution of all simulations equally. Nevertheless, we 
find most optimistic conditions for the onset of an explo- 
sion in any case, in particular for the 15 M , using LS180 
and LS375. Whether this remains true for LS including 
these yet missing contributions, will have to be shown in 
a future exploration. On the other hand, such optimistic 
conditions were not obtained in ^600 ms postbounce for 
our 2D simulations using SHEN, even without the above 
mentioned important sources of energy loss. It indicates 
that it may be much more difficult to obtain neutrino- 
driven explosions for this progenitor model using SHEN 
EOS. 

4.3. Results of a model with 11.2 Mq 

In addition to the 15 Mq progenitor from 
iWooslev fe Weaver! ()1995|) . we perfor med 2D simu- 
lation s of an 11.2 Mq progenitor from IWooslev et al.l 
(|2002f) . The 11.2 Mq star has been used in several 
studies before, w here neutrino-driven explosions were 
obtai ned in 2D dBuras et all l2006al iMarek fc Jankal 
2009) and 3D (jTakiwaki et al.ll2012t ) simulations. Figure 



20(a)| shows the time evolution of the average shock 
radius in our 2D models, comparing LS180 (red solid 
line) and SHEN (grey dot-dashed line). For the 15 
Mq progenitor model discussed above, only a passive 
shock expansion was obtained due to the decreasing 
ram pressure of LS EOS. In contrast to this, both EOS 
lead to the successful launch of an explosion in the case 
of the 11.2 Mq progenitor, for which the expansion 
of the shock is acc ompani ed by mass outflow. This is 
illustrated in Figure 20(b)[ which shows the post-bounce 
evolution of selected mass shells for the simulation 
using SHEN. This Figure should be compared with 
Figure [SI which shows the 15 Mq model using LS180. 
The differences are caused by the different progenitor 
structures of the 15 Mq and the 11.2 Mq progenitor 
models. The latter has a much steeper density gradient 



at the interface between the iron core and the Si-layer, 
as well as between the Si-layer and the C/O-layer. A 
more detailed analysis regarding the criteria for the 
successful launch of a neutrino-d ri ven ex plosion has 
already been given in iBuras et all (|2006af ). where the 
same two progenitor stars have been compared. 

For both, SHEN and LS180, the post-bounce accre- 
tion phase lasts for only about 50 ms. After this short 
accretion phase, the shock waves start to expand contin- 
uously to larger radii. This is significantly earlier tha n 
in the simul a tions r eported by Marek fc Jankal (|2009| ): 
IMiiller et all (|2012al ). In the latter, the explosion was 
delayed by several 100 ms. Our models are not taking 
general relativistic effects into account and neglect the 
emission of fi— and r-ncutrinos. Both effects contribute 
toward a larger radius of the shock in the post-bounce 
phase compared to models with standard input physics. 
Moreover, the findings described above in the compari- 
son between SHEN and LS models of the 15 Mq progen- 
itor are seen in the 11.2 Mq model as well. The models 
based on SHEN lead to less favorable conditions for the 
explosion. The difference between SHEN and LS180 in 
the 11.2 Mq model reduces to a slight delay in the on- 
set time of the explosion and a slower shock expansion. 
As discussed above, it is due to less efficient neutrino 
heating in the SHEN case. Nevertheless, both explosions 
are weak. At the end of the simulations, the diagnostic 
energy is i^dia ~ 10 50 erg. 

A snapshot of the Y e and entropy per baryon distri- 
butions is shown in Figure [21] at 555 ms post bounce. 
At this this is the last candidate, next esc will revert 
to uncompleted text, ime the model has already entered 
the neutrino-driven wind phase, in which regions of con- 
stant entropy per baryon emerge around the protoneu- 
tron star (see fig. 121"]) . Based on spherically symmetric 
simulations o f massi v e iron -core collapse and explosions, 
iFischer et all (120101 120111) c o nfirm ed recently the ex- 
pectation of lOian fc Wooslevl (|1996f ). that the neutrino- 
driven wind attains proton-rich conditions. Our models 
are consistent with these results: During the early wind 
phase we obtain Y e ~ 0.5 — 0.55 in the ejected matter. 

5. SUMMARY AND DISCUSSION 

We performed ID and 2D radiation hydrodynamic sim- 
ulations that include spectral neutrino-radiation trans- 
fer. The simulations are launched from a 15 Mq pro- 
genitor star using four different EOS, corresponding to 
three different incompressibilities of the LS EOS and one 
SHEN EOS. LS is based on a non-relativistic approach 
and has a symmetry energy of 29.3 MeV. The incom- 
prcssibility parameter can be set to 180 (LS180), 220 
(LS220), or 375 MeV (LS375). SHEN is based on rela- 
tivistic mean field theory and the Thomas-Fermi approx- 
imation. It has a symmetry energy of 36.9 MeV and an 
incompressibility of 281 MeV. This set of of EOS repre- 
sents a wide range of nuclear matter properties. 

In ID, none of the simulations produced an explo- 
sion for the considered simulation times up to 1 s post 
bounce. The observed similarities between LS180, LS220 
and LS375, as well as differences to SHEN, concerning 
the conditions at bounce, the propagation of the shock, 
and the evolution of neutrino luminosities and mean en- 
ergies, are in agreement with previous studies. The 
incompressibility and symmetry energy of the nuclear 
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Figure 19. Post-bounce evolution of the heating efficiency (left) and the growth parameter of convection (right). In the right panel, the 
dotted horizon line represents the critical value of 3. Both quantities do not exhibit the significant EOS dependence. 



3 

M 

o 
o 

C/5 



2500 



2000 



1500 



1000 



500 



:' us 180 — — 




/ - 


Shen 




/ - 






/ ^ 

/ / : 






/ / 






/ / - 






/' ; 




//" 





-100 100 200 300 

Time after Bounce [ms] 

(a) Average shock radius. 



400 




150 200 250 300 350 
Time after boune [msl 
(b) The mass trajectories. 



500 



Figure 20. Post-bounce evolution of the average shock radius (left) and selected mass elements (right) for a 2D simulation of a 11.2 Mq 
star. In the right panel, the thick black line is the angle-averaged shock radius. The thin black and grey lines show mass elements with an 
enclosed mass from f.O to 1.4 Mq (thin black) and from 1.3 to 1.4 Mq (grey) at intervals of O.Of Mq (thin black) and O.OOf Mq (grey), 
respectively. 



EOS depend on the density and determine the condi- 
tions from which the ex plosion may emerge. However, 
as recently discussed in iSteiner et al.l (|2012f ). it is not 
possible to reduce the characteristics of the EOS to a 
single important parameter. Note that the re are several 
constraints on the nuclear p roperties (e.g.. ISteiner et all 
l2010tlLattimer k Liirj|2012j) . These suggest that, among 
the employed EOS in this study, SHEN has too high 
symmetry energy, LS180 has too low incompressibility, 
and LS375 has too high incompressibility. Only LS220 is 
compatible with terrestrial experiments and the 1.97 Mq 
Demorest et al. pulsar mass constraint for cold neutron 
star matter. 

In 2D, the differences between models with LS and 
SHEN EOS are consistent with the differences between 
the ID models. For example, we find significantly more 
favorable conditions for the launch of a neutrino-driven 
explosion in the models using the LS EOS: The shock ex- 
pands continuously to large radii, which are not reached 
in the models using SHEN at comparable times. We an- 



alyzed these different post-bounce evolutions and found 
that models with LS show a highly turbulent velocity 
field producing prompt convection during the early post- 
bounce phase. Moreover, LS leads to to higher electron 
antineutrino luminosity and more efficient neutrino heat- 
ing. More mass is accumulated in the gain region, which 
leads to a higher entropy per baryon at the early post- 
bounce phase. In addition, the standing accretion shock 
instability (SASI) pushes the standing accretion shock 
to increasingly larger radii after a post-bounce time on 
the order of 100 ms. In contrast to these models based 
on LS, the models based on SHEN show weak neutrino 
heating and the standing accretion shock is only oscil- 
lating around a mean radius for simulation times up to 
600 ms post bounce. Note, however, that our simula- 
tions ignore inelastic scattering on electrons, which is 
important during the collapse phase, and the emission 
of heavy-lepton neutrinos, which contributes to the en- 
ergy loss during the post-bounce phase. These reactions 
would both make our models more pessimistic, so that 
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Figure 21. The Y e (left panel) and entropy per baryon (right 
panel) distributions of a 2D simulation using SHEN. Shown is a 
snapshot at the post-bounce time oft 555 ms, when the model has 
already entered the neutrino-driven wind phase. 

the 2D models based on SHEN appear unlikely to pro- 
duce an explosion with the complete set of neutrino reac- 
tions. Hence, the different properties of the nuclear EOS 
can explain some apparently contradictory results found 
in previous studies, for example the outcome of neutrino - 
driven explos i ons in s imulations of iBuras et al.l (|2006aD ; 
IBruenn et all (12001) : ISuwa et al.1 ([20100 using LS EOS 
and the absence of neutrino- driven explosions in simula- 
tions of lBurrows "et^lpOOfl using SHEN. Very recently, 
ICouchl (|2012l ) investigated the EOS dependence based 
on simulati ons with parameteri z ed neu ti rno luminosities 
simila r to iMurphv fc Burrows! (|2008l ); iNordhaus et al.l 
(|2010fklHanke et al.l (|2012D . He found that a 15 M Q star 
is more difficult to explode using SHEN than LS EOS, 
which is also consistent with the result of this paper. 

The most optimistic model, with respect to a possible 
neutrino-driven explosion, was obtained using LS180, for 
which we explore different indicators for the possible on- 
set of an explosion. We note again that the neutrino- 
driven shoc k revival occurres much faster in this model 
compared to lMarek fc Ja nka (2009|) because of the afore- 
mentioned neglect of the additional cooling processes. 
The energy for the shock expansion stems from aspheri- 
cal mass accretion from large radii deep onto the surface 
of the protoneutron star, generating a semi-stationary 
advection cycle. Thus, matter which expands behind 
the shock wave experiences a continuous flux of neutri- 
nos from the accreted matter. It remains to be shown 
whether the asymmetry of the accretion pattern is simi- 
larly pronounced in three-dimensional supernova simula- 
tions. A small difference between the v e and v e spectra 
leads to proton-rich conditions. Whether these findings 
will remain after having included a treatment of charged- 
current processes th at is consistent with the EO S, as 
discussed recently in iMartmez-Pinedo et al.l (|2012t ) and 
IRoberts fc Reddvl (12012D . remains to be shown in fu- 
ture studies. Such improvements may become important 
especially at later post-bounce times that we have not 
reached in the present simulations. 



Based on the described explosion model we evaluate 
different indicators that have been suggested in the lit- 
erature to diagnose the onset of an explosion during the 
post-bounce phase. We confir m that mass- weighted aver - 
age entropy in the gain region (jMurphv fc Burrows|l2008D 
and the mass enclosed in the gain region itself (jJankal 
|2001[ ) both show the expected difference between the op- 
timistic LS and the pessimistic SH EN models. With re- 
spect to the ante-sonic condition (|Peicha fc Thompson! 
I2012T ). we find that it applies as well. However, fast con- 
vection can lead to transient peaks in the indicator that 
may not be interpreted as a possible onset of the explo- 
sion. Also the ratio between advection time scale and 
heating time scale has been su ggested as explosion in- 
dicator dThompson et al.l 120051) . This indicator is also 
affected by multidimensional effects and one would have 
to carefully specify how the time scales have to be av- 
eraged in a multidimensional setting. The evaluation of 
the time scales based on angularly integrated quantities 
does not lead to a robust indicator for the onset of an 
explosion. Also the evolution of the radius of the pro- 
toneutron star does not serve as a reliable predictor of 
the explosion. 

In order to extend our EOS comparison study to a 
broader range of initial progenitor models, and in par- 
ticular to confirm the optimistic results obtained using 
LS, we also include a lower mass iron-core progenitor of 
11.2 Mq into our investigation. For this model, we ob- 
tain neutrino-driven explosions for both LS and SHEN 
EOS. The differences observed between both simulations 
remain qualitatively similar to the ones discussed above 
for the 15 Mq models. In comparison to LS, SHEN leads 
to a slightly delayed onset of the explosion with a smaller 
explosion energy. 

However, it should be noted that the simulations in 
this paper are only a very first step towards more realis- 
tic supernova models. In addition to the simplifications 
of weak processes and the omission of heavy-lepton 
neutrinos, general rclativistic effects should be consid- 
ered. Additionally, the ray-by-ray approximation may 
lead to an overcstimation of the directional dependence 
of neutrino anisotropics and of radial oscillations of 
SASI. A full-angl e transpor t will give us a mo re re- 
fined answer (see lOtt et all 120081: iBrandt et all 120111: 
iSumivoshi fc Yamadall2012l l Moreover, due to the coor- 
dinate symmetry axis, the SASI develops preferentially 
along this axis. It could thus provide more favorable con- 
dition for the growth of I = 1 mode of the SASI and also 
for the possible onset of explosions along this direction. 
In the appendix we point to the clear and crucial impact 
the angular width of the computational domain has on 
the activity of the SASI, and in consequence on the 
presence or absence of the explosion. Hence, supernova 
modelers are gra dually switching gears from 2D to 3D 
simulations ( e.g.. Ilwakami et al| |2008: Sc heidegger et all 



2008t iKotakeet al.l 120091: iScheidegger et al 



Nordhaus et all 20101: Wongwathanarat et al 



Takiwaki et alJ I2012t lHanke et al.l 120121: iBurrows et all 



20121 ) . in which the SASI and convection were observed 



to develop much more stochastically in 3D than in 2D. 
This shows that complementary 3D supernova models 
are indeed necessary to pin down the impact of the 
EOS on the neutrino-driven mechanism. Together with 
the 3D effects, it would also be interesting to study 
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Figure 22. Post-bounce evolution of the average shock radius for different mesh configuration. 



the possible EOS impacts on the gravitational- wave 



emission (e.g.. 


Marek et al.H2009L 


iKotake et al.l 


20121 for a recent 



review). We will address these questions in forthcoming 
papers. 
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APPENDIX 
IMPORTANCE OF SASI 

In addition to the full it simulation, we perform simulations with 0° < 8 < 90° and 45° ^ 8 < 1 3 5° for LS180 , whic h 
shows SASI activity . The latter grid setting was often used in the past, cf., iBurrows et al.l (|1995l ); iFrver et ail (11999ft; 
Buras et al.1 (|2006al ). The SASI is thought to be fastest growing for low-£ modes, such as i =1 and 2 (jBlondin et al.1 
20031 ). The restriction to 0° < 8 < 90° does not include I = 1, but it includes the i = 2 mode, while the restriction 
to 45° < 8 < 135° does include none of them. Figure [22] shows the post-bounce evolution of the shock average 
radius for different simulations using LS180 and the three above mentioned resolutions. The red-solid, green-dotted, 
and blue-dashed lines indicate the full ir simulation, tt/2 for the northern pole, and 7r/2 around the equatorial plane, 
respectively. We find that the simulations with 0° < 8 < 90° represent an optimistic case with respect to possible 
explosions with continuous shock expansion to larger radii. For 45° < 8 < 135°, the shock wave does not continue to 
expand to larger radii and represents the least optimistic scenario with respect to possible explosions. For the latter 
case, where the development of possible SASI activity is suppressed due to the chosen angular resolution, possible 
explosions cannot be expected. Thus, large degrees of freedom with respect to angular resolution, allowing for the 
evolution of SASI, is also an important requirement being able to make predictions about the possible onset of an 
explosion. 
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